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Abstract 

In models with gauge mediated supersymmetry breaking the gravitino is the 
lightest supersymmetric particle. If R-parity is violated the gravitino decays, but 
with a half-live far exceeding the age of the universe and thus is, in principle, 
a candidate for the dark matter. We consider the decays of the next-to-lightest 
supersymmetric particle, assumed to be the neutralino. We show that in models 
where the breaking of R-parity is bilinear, the condition that R-parity violation 
explains correctly the measured neutrino masses hxes the branching ratio of the 
decay Xi Gj in the range 10“^ — 10“^, if the gravitino mass is in the range 
required to solve the dark matter problem, i.e. of the order (few) 100 eV. This 
scenario is therefore directly testable at the next generation of colliders. 
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1 Introduction 


Models of gauge mediated supersymmetry breaking (GMSB) generically pre¬ 
dict that the lightest supersymmetric particle (LSP) is the gravitino [1]. Such a 
light gravitino, in principle, is a candidate for the non-baryonic dark matter of 
the universe [2], The smallness of the gravitino couplings, however, make such 
a scenario extremely difficult to test. Direct detection of gravitino dark mat¬ 
ter in scattering experiments or indirectly via decays/annihilation to gamma 
rays is hopeless [3] and consequently gravitino dark matter has received rather 
scarce attention. 


The purpose of the present letter is to show that in models with bilinear 
breaking of R-parity the branching ratio of the decay of the neutralino into a 
gravitino and a photon is fixed by data on neutrino masses up to a factor of 
771^^2) if ffi® neutralino is the next-to lightest supersymmetric particle (NLSP). 
A measurement of this branching ratio thus implies a “measurement” of the 
gravitino mass, 7773 / 2 . Approximate knowledge of 7773/2 in turn can be used to 
constrain the conjecture that the gravitino is the (major component of the) 
dark matter in the universe. 


On the one hand, supersymmetric models with R-parity violation can ex¬ 
plain [4-7] current data on neutrino masses and mixings [8] without invoking 
any GUT-scale physics. On the other hand, in supersymmetric models with 
R-parity violation the LSP decays. For all superpartners of standard model 
particles these decays proceed at rates that even the most tiny amount of R- 
parity violation rules out MSSM particles as dark matter. A light gravitino, 
however, couples so weakly to standard model particles that its half-live far 
exceeds the age of the universe even for R-parity violating couplings as large 
as 0{1) [9,10], see also section 2. Thus, contrary to popular believe, supersym¬ 
metry holds the promise to solve the dark matter problem even if R-parity is 
violated. 


If gravitinos were in thermal equilibrium in the early universe (and assuming 
that there is no non-standard physics between gravitino decoupling and the 
time of nucleosynthesis, see below), the contribution of gravitinos to the matter 
content of the universe can be estimated to be [2] 


123/2/7^ 





^3/2 

100 


( 1 ) 


Here, 123/2 is the density of gravitinos in units of the critical density, h is the 
Hubble parameter in units of 100 km s“^ Mpc“^ and is the effective number 
of degrees of freedom at the time of gravitino decoupling. Depending on the 
so-far unknown supersymmetric particle spectrum one expects gf* ~ 90 — 140 
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[11]. Current data give the matter density of the universe [12] as VLuh? — 
0.134 ± 0.006, from which ~ 0.023 ± 0.001 is in the form of baryons. 

Particle dark matter (DM) is usually classihed according to its free-streaming 
length [13] as either “hot”, “warm” or “cold” DM. There is a general con¬ 
sensus that hot DM is ruled out [3,12,13]. Cold DM is usually considered the 
best choice [3,13] to £t large-scale structure data. However, on galactic and 
sub-galactic scales pure cold DM seems to produce too much power, see for 
example [14,15] and references therein. ^ To resolve the deficits of cold DM, 
some groups considered warm DM variants, claiming that WDM does actu¬ 
ally provide a better fit [15,17]. However, constraints on the free streaming 
length (and thus the mass) of WDM particles can be derived from data of the 
Lyman-o; forest [18], and a lower limit of mwoM > 0.55 keV for thermal relics 
is quoted in [19]. 

The lower limit on the mass of WDM particles given in [19] seems to be in 
conflict with the conjecture that gravitinos are the dark matter, see Eq.(l). 
However, in the derivation of the gravitino density [2] it is assumed that the 
universe has a “standard” thermal history. Producing additional entropy af¬ 
ter the time of gravitino decoupling would dilute the density of gravitinos, ^ 
compared to the estimate Eq.(l), two variations of this idea are discussed in 
[20,21]. Both [20] and [21] assume that entropy is produced by the “late” de¬ 
cay of messenger particles. ® Baltz and Murayama [20] argue that the lightest 
messenger particle might decay through an intermediate heavy particle with 
mass mx — 10^^ GeV, which leads to a messenger decay width sufficiently 
small to dilute the gravitino density by a factor of (5 — 8). Fujii and Yanagida, 
on the other hand, claim that adding a constant messenger number violat¬ 
ing term to the superpotential, messenger widths of the “correct” order of 
magnitude are naturally obtained [21]. 

Given this discussion, we think it is fair to say that gravitinos with a mass 
in the range of (P(O.l)-(P(l) keV are interesting dark matter candidates. Gon- 
straining the gravitino mass to be much smaller than given by Eq.(l) would 
rule out gravitinos as DM. 

NLSP decays in R-parity violating variants of GMSB have been considered 
previously [24]. For the case of bilinear R-parity breaking the authors of [24] 
point out that a bound on |A|/^/detAf^o of (very roughly) the order 0(10“®) 

^ For a different point of view, see for example [16]. 

® The other logical possibility, i.e. to raise to values of the order of ~ (600 — 
700) by the introduction of a sea of new particles, does not seem very economical. 

® It has been speculated that messenger particles themselves might provide the 
dark matter, see for example [22,23]. However, messengers tend to overdose the 
universe unless the mass of the lightest messenger is rather low, of the order of 
miMP - 0(1) TeV. 
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for ^/F = 10® GeV can be obtained from the requirement r(x? —^ Gj) > 
X]r(x? I/ip). Fits [4] to current neutrino data [8] require similar, although 
somewhat larger, values for bilinear R-parity violating parameters, see next 
section. In our numerical calculation we thus hnd r(x® — G7)/X]r(x5 — Ifip 
) < 1, unless the gravitino mass is much smaller than indicated in eq. (1). 

The remainder of this paper is organized as follows. In the next section we 
present some approximate formulas for the decay of the neutralino NLSP. 
These estimates serve to understand the results of the numerical analysis, 
presented next. We then close with a short summary. 


2 Semi-analytical estimates 


In this section we give some (semi-) analytical formulas for the decay of a neu¬ 
tralino NLSP. This will facilitate the understanding of our numerical results 
presented below. In GMSB with R-parity violation, a neutralino NLSP can ei¬ 
ther decay into a gravitino and a photon or via R-parity violating interactions 
directly to standard model particles. 

We take into account only bilinear R-parity violating (BRpV) terms, namely. 


W — Wmssm + ^iLiHu, 

Rsoft = (2) 

Eq. (2) can be considered as a minimal model of R-parity violation. The new 
terms in Rjoft induce vacuum expectation values for the scalar neutrinos Vi. 
One can either treat R* or Vi as free parameters of the model, since they are 
connected by the tadpole equations. 

In models with bilinear breaking of R-parity, one neutrino mass is generated 
at tree-level, while the other neutrino masses are due to 1-loop corrections 
[4]. The decay width of the NLSP to standard model particles is related to 
the neutrino masses, thus we first discuss some approximate formulas for the 
calculation of neutrino masses in BRpV models. 


The tree-level contribution of BRpV to neutrino masses is given as 


m 


tree 


rriF^ 

4detAI^o 



(3) 


Here, is the “photino mass” = g'^Mi + g'‘^M 2 , detAl^o is the deter¬ 
minant of the (4 X 4) MSSM neutralino mass matrix and A is the so-called 
alignment vector. A* = eiVd + Uj/i. 
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The dominant 1-loop corrections to the neutrino mass matrix are usually due 
to bottom/sbottom and tau/stau loops and are, very roughly, of order [4] 






(3hlsm{2ei)mbABy^^ + hi sm{29r)m-, 




(g? + gj) 

/i2 


( 4 ) 


Here, hb {hr) are the bottom (r) Yukawa coupling, 9^ {9f) is the mixing angle 
in the sbottom (stau) sector, ABq^ is the difference of two Passarino-Veltman 
Ho"functions, essentially AB^^ ~ ln(m^/m^) for sfermion masses much larger 
than the corresponding fermion masses. And, hnally e are the superpotential 
parameters e rotated to the basis where the tree-level neutrino mass matrix is 
diagonal. 

Eqs (3) and (4) produce a hierarchical neutrino spectrum. We will assume 
that the tree-level contribution is larger than the 1-loop correction. Thus, we 
identify ~ ^Am%^ ~ 0.04—0.06 eV and p ~ Am q ~ 0.009 eV. For 
any given choice of R-parity conserving SUSY parameters then |A|/y^dutA^ 
and \e\/fi are approximately hxed by neutrino masses. Typical values are [4]: 
|A|/ydetAf^o ~ (few) 10“® and |e|//r~ (few) 10“"^. 

The neutralino will decay to three SM fermion final states or, if kinematically 
allowed, into gauge bosons and leptons, and To estimate the 

most important decay widths, we will make use of the approximate neutralino 
couplings in first order expansion in small parameters as given in [25]. 

Consider Erst the decay to gauge bosons. In GMSB scenarios the lightest 
neutralino is usually bino dominated [1]. Binos couple to gauge bosons pro¬ 
portional to Aj. With couplings from [25] and Eq.(3) we estimate 


i 


g^g'^M2m^o 

{IQuMimrl) 


/(m^/m|o)m' 


Tree 

V • 


( 5 ) 


Here, f{x) is a phase space factor, given hj f{x) = ^ ~ ^ A similar 

expression holds for r(x5 — Z^^i) with an additional prefactor of l/(4c^). 
Assuming Mi ~ M 2/2 and \fj,\/Mi ~ 4, as is typical for GMSB models, results 
very roughly in T ~ 2 x 10’ eV. 

Neutralino decays to three fermions can also be mediated by scalar quark and 
scalar lepton exchange. With approximate scalar lepton couplings from [26] 


^ The decay to a light Higgs plus neutrinos, Pv, is also possible. However, for 
generic GMSB parameter choices it is less important. 
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we estimate that the decay r(x 5 ’ —^ ur^r^) is very roughly of order 


n -L /^\2/Tfl^0\A Tft~o 

r(x. - P) ~ ^(-) (^) 9(^)mi; 


m^o - 
rrif 


m- 

A 

m% 


( 6 ) 


Here, 

g{y) = ^(-5/2 + ?,/y + (-1 + l/y){-l + 3/|/)ln(l - y)) (7) 

y 

and we have normalized g{y) conveniently such that 5 '( 0 ) —> 1 , thus g{y) varies 
between [1,6]. Note, that r(x'j’ —goes to zero proportional to the 
fourth power of m^o/rrif. Since in GMSB the (right) scalar tan is never very 
much heavier than the neutralino, contrary to squarks and other scalar leptons, 
Eq.( 6 ) usually dominates over other Feynman graphs with scalar exchange. 

Eq.( 6 ) can take a wide range of values. Just to give a flavour of its typical 
size, for ~ 3 X 10“^, tan/? = 10, and m^o = 100 GeV and = 120 GeV 
one hnds r(x5 —> ~ 3 x 10“^ eV. We have checked numerically, see 

next section, that the decays described by eqs (5) and ( 6 ) are usually the most 
important RpV decay channels in GMSB. 

The decay width of a neutralino NLSP to gravitino-photon is given by [1] 


r(x? - Gy) = 






-1.2 X 10“®fi;2( 


/ m^o xo/ iuuev x 

Vino GpV/ V ) 


5/100 eV\2 


100 GeV 


eV 


m3/2 


( 8 ) 


Here, ® = | cos OwNu + sin 6 *vi/iVi 2 1. Neutralinos can also decay into Xi 

GZ^ and y? — Gh^ and we include those channels in our numerical calcu¬ 
lation. However, these hnal states are usually less important than Xi Gj 
and also do not give a promising signal. We will therefore not discuss them in 
further details. 


From eqs (5), (6) and (8) we can very roughly estimate a branching ratio of 
Br{xi Gy) — 10“*^^“^^ for 7713/2 of G(IOO) eV. This is the main result of 
the current paper. We will back up this estimate with a numerically exact 
calculation in the next section. 


We note for comparison that if R-parity is conserved, Eq.(8) gives a typical 
decay length of cr — 20(^^^)^ m, for a neutralino mass of m^o — 100 GeV. 

® In MSSM notation. In the notation of [4]: Nij -p N^j 
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According to [27] the ATLAS detector at LHC should be able to measure even 
such a “large” decay length rather well given sufficient luminosity ( 100 / 6 “^). 


Finally, we have checked that the gravitino itself lives long enough to be a dark 
matter candidate. Following [10] the decay width of G —z/y can be calculated 
from the photino content of the neutrino as 


i 


1 


\U^u 


2 ^ 1/2 

Mpi 


(9) 


Here, = Z]i=i I cosOwNn + sm9wNi2\‘^. The coupling matrices Nn and 

Ni 2 can be calculated perturbatively [4] and are approximately hxed from the 
neutrino masses. For Mi = 100 GeV and 7713/2 = 100 eV, we hud ~ 

3.5 X This corresponds to a half-life of the order of 10^^ Hubble 

times. 


3 Numerical results 


We stress that none of the approximations discussed above were used in our nu¬ 
merical analysis. Numerical results presented below have been obtained as fol¬ 
lows. We generated supersymmetric particle spectra using the package SPheno 
[28]. GMSB is characterized by two mass scales, the scale of supersymmetry 
breaking F and the messenger mass Mm- Since F is related to the gravitino 
mass via [ 1 ] 


k-/3Mp 


( 10 ) 


we trade F for the gravitino mass and vary = F/Mm independently. 

For details see [1]. For dehniteness we take k = 1/20. We have checked that 
our results depend only very weakly on the exact value of k as long as /c < 1. In 
addition, we have tan (3, the sign of fi and the number of messengers, 775 , as free 
parameters. We use only 775 = 1,2 because for larger values of 775 the neutralino 
is rarely the NLSP, since scalar masses scale as nif/m^o ~ l/^/n^ce — b.ce-A 

We check the generated spectra for a number of phenomenological limits [12]: 
777 -+ > 105 GeV, mr, > 95 GeV and nif > 82 GeV. Lower limits on the mass of 
the lightest neutralino in GMSB have been published by all LEP collaborations 
[29], bounds are between 92—100 GeV depending on the details of the analysis. 
However, we have found that the most important constraint for us is the lower 
limit on the mass of the lightest Higgs. [30], which essentially cuts out all points 
with m^o < 100 GeV from our scans. 
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R-parity violation is then included into SPheno [28] as described for neu¬ 
trino masses in [4] and for neutralino decays in [25]. Special care is taken to 
diagonalize the neutrino-neutralino mass matrix at 1 -loop order to give neu¬ 
trino masses and mixings compatible with the values indicated by atmospheric 
and solar neutrino experiments [ 8 ]. For each set of GMSB parameters, m 3 / 2 , 
fyGMSB^ tan (3, sgn{fj,) and ns, this results in a restricted range of A* and e*, as 
discussed above, from which then the RpV neutralino decays are calculated. 
From the discussion in the previous section one expects that the errors in our 
calculated branching ratios scale proportional to the errors in the neutrino 
masses and we have checked numerically that this assertion is correct. Results 
discussed below use neutrino masses near the best £t points for solar and at¬ 
mospheric physics [ 8 ]. Finally, for the calculation of the branching ratio into 
gravitino plus photon Eq. ( 8 ) is used. 

Our numerical results show that the dominant hnal states are usually either 
or (and Z^u). Gauge boson hnal states become more important 

the larger rrif — m^o, as discussed above, and for larger m^o. We hud that for 
> 150 GeV the hnal state Pu can reach a branching ratio of the order 
of 5 — 15 %. All other hnal states usually have branching ratios which are 
smaller. Especially, we hnd that the hnal state bbu is less important than in 
an mSugra scenario [25]. This can be understood as being due to the smaller 
ratio of in GMSB compared to mSugra. With typical total widths 

in the range of roughly (9(10“*^) — 0(10“^) eV, we expect that the neutralino 
decays with a displaced vertex. 

In Figure 1 we show the branching ratio BR{'^ G'y) as a function of NLSP 
mass for ns = 1, two values of tan/3, both signs of /i and for a fixed value of 
'nT' 3/2 = 0.1 keV. The branching ratio rises strongly with increasing neutralino 
mass, as expected from eq. ( 8 ). The plot also shows that the dependence on 
tan (3 is rather weak, changing tan (3 from 10 to 35 induces a change in the 
branching ratio up to a factor of ~ 2. It is also obvious that the sign of p is 
not decisive. Ghoosing ns = 2 reduces the branching ratio by typically a factor 
of ~ 3 compared to the results shown. 

The experimental signal for the hnal state Gy is 7 . In R-parity violating 
models the neutralino has another decay mode which gives the same exper¬ 
imental signal, namely ^7- BRpV this occurs at 1-loop order. To 

estimate this background we have done a calculation of > z/ 7 ). Fig. 

( 2 ) shows BR{x\ G'y)/BR{xi z^y) as a function of the gravitino mass, 
for two diherent choices of tan (3 and for ns = 1 for various diherent values of 
m^o. The ratio depends strongly on the neutralino and gravitino masses, for 
^ 3/2 < 500 eV it is always larger than 1. For neutralino masses greater than 
about m^o ~ 150 GeV z/y never seems to be a serious problem up to gravitino 
masses 1713/2 = 2 keV. 



T 

O i-H 

05 



Fig. 1. BR{x^ —> G 7 ) as function of the lightest neutralino mass, m^o [GeV]. 
Full lines are for /i > 0, dashed lines fi < 0. Light (on colour printers magenta): 
tan/3 = 10, Dark (blue): tan/3 = 35. 



Fig. 2. Ratio BR{xi —> G'y)/BR{xi —> as function of m ^/2 [eV]. Full lines are 

for tan/3 = 10, dashed lines for tan/3 = 35. The plots shows the case ns = 1, for 
ns = 2 the ratio is typically a factor (2-3) smaller. The different lines are for (from 
bottom to top) m^o = 100 — 500 GeV in steps of 100 GeV. 
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Finally, fig. (3) shows our main result. Here we plot BR{x^ Gy) as a 
function of the gravitino mass for different values of m^o, two values of tan (3 
and ^5 = 1 . 77.5 = 2 leads to branching ratios approximately a factor of up to 
3 smaller. Depending on the neutralino mass, > Gy) is larger than 

10“^ for values of 7713/2 — 0.5 keV up to 777 . 3/2 ~ 2 keV. At the LHC one 
expects to produce very roughly of the order of G(10®) - G(IO^) events from 
supersymmetry, depending mainly on squark and gluino masses. Thus we think 
that sufficient statistics to measure branching ratios as small as 10“^ should be 
possible. We conclude therefore, that for cosmologically interesting ranges for 
the gravitino mass measurably large branching ratios BR{xi —>■ Gy) should 
exist. 



Fig. 3. Ratio BR{xi —> Gy) as function of the gravitino mass, 7773/2 [eV]. Full lines 
are for tan/3 = 10, dashed lines tan/3 = 35. Lines from bottom to top are for 
different values of rn^o, m^o = 100 — 500 GeV in steps of 100 GeV. BR{x\ —> Gy) 
is larger than 10“^ for values of 7773/2 between 0.5 — 2.0 keV, depending on the 
neutralino mass. 
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4 Conclusions 


We have discussed decay properties of the lightest neutralino in models with 
gauge mediated supersymmetry breaking in which the neutrino masses and 
mixings are explained by bilinear R-parity violation. Once the BRpV param¬ 
eters are approximately hxed with information from the neutrino sector, the 
branching ratio into gravitino plus photon is hxed to be in the range 
for a gravitino mass of the order of (few) 100 eV. The branching ratio decreases 
with increasing gravitino mass. In the scenario discussed one can therefore test 
whether the gravitino gives a signihcant contribution to the dark matter of 
the universe by a “simple” counting experiment. 

Can one do better? - Concerning the gravitino mass, the answer is yes. A 
measurement of the decay length of the neutralino would hx the total neu¬ 
tralino decay width, independent of our assumptions about neutrino masses. 
Knowing the width and the mass of the neutralino hxes the gravitino mass 
from the measurement of G'j) in a much tighter range than what we 

have been able to do. However, one has to admit that even knowing the grav¬ 
itino mass rather well, one can not calculate ^ 3/2 reliably from 7713/2 without 
making specihc assumptions on the thermal history of the universe. In this 
sense, gravitino DM can be ruled out by the measurement we have discussed, 
but never “experimentally conhrmed”. 
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